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^ ■ ABSTRACT 
^ . 

CN ! Of the cosmological gamma-ray bursts, GRB 011121 has the lowest redshift, 
z = 0.36. More importantly, the multi-color excess in the afterglow detected 

^ ! in the Hubble Space Telescope (HST) light curves is compelling observational 

\^ i evidence for an underlying supernova. Here we present near-infrared and radio 

^ i observations of the afterglow. We undertake a comprehensive modeling of these 

O ! observations and those reported in the literature and find good evidence favoring 

^ I a wind-fed circumburst medium. In detail, we infer the progenitor had a mass 

"rH i loss rate of M ~ 10~''/f^3 Mq yr~^ where f^s is the speed of the wind from the 
progenitor in units of 10^ km s ^. This mass loss rate is similar to that inferred 
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for the progenitor of SN 1998bw which has been associated with GRB 980425. 
Our data, taken in conjunction with the HST results of Bloom et al. (2002), 
provide a consistent picture: the long duration GRB 011121 had a massive star 
progenitor which exploded as a supernova at about the same time as the GRB 
event. 

Subject headings: gamma rays: bursts 

1. Introduction 

On 2001 November 21 at 18:47:21 UT, GRB 011121 was detected and localized by the 
Italian-Dutch satellite BeppoSAX (Piro 2001). The localization was further improved by 
the Interplanetary Network (Hurley et al. 2001) and an optical transient was identified by 
the OGLE group (Wyrzykowski, Stanek & Garnavich 2001). Spectroscopy of the transient 
revealed emission lines interpreted as arising from the host galaxy at a redshift oi z — 0.36 
(Infante et al. 2001). 

Low redshift GRBs are particularly valuable in uncovering the origin of GRBs. If GRBs 
result from the death of massive stars then it is reasonable to expect an underlying supernova 
(SN). Bloom et al. (1999) attributed a late-time red excess seen in the afterglow emission 
of GRB 980326 to an underlying SN. This result triggered searches for similar excesses with 
no clear success save GRB 970228 (Reichart 1999; Galama et al. 2000). The low redshift is 
critical to such searches since the SN contribution is expected to exhibit strong absorption 
below 4000 A (see Bloom et al. 1999). 

Given this motivation, we triggered a sequence of multi-color and multi-epoch Wide 
Field and Planetary Camera 2 (WFPC-2) observations with the Hubble Space Telescope 
(HST). Garnavich et al. (2002) noted that the i?-band flux of the first epoch of the HST 
observations was significantly in excess of the extrapolation of the power law decay of the 
early ground-based optical afterglow and attributed this to an underlying SN component. 
In Paper I (Bloom et al. 2002) we presented four-epoch multi-color HST light curves and 
show the data are explained by an underlying SN similar to SN 1998bw (Galama et al. 
1998) except fainter by about 2/3 magnitude. At this point, there appears to be compeUing 
evidence for GRB 011121 to be associated with a SN which exploded at about the same time 
as the gamma-ray event (Bloom et al. 2002). 

This GRB-SN hnk is an essential expectation in the coUapsar model (Woosley 1993) in 
which GRBs result from the death of certain massive stars. Another essential consequence of 
any massive star origin for GRBs, as noted by Chevalier & Li (1999), is a circumburst medium 



-3- 



fed by the inevitable and copious mass loss suffered by massive stars throughout their lives. 
Afterglow observations are excellently suited to determining not only the geometry of the 
explosion but also the distribution of circumburst matter. Unfortunately, until now there 
has been no clear evidence for a wind-fed circumburst medium (density, p oc with s ~ 2; 
here r is the distance from the explosion site) in the afterglow of any cosmologically located 
GRB. 

Here we report near-infrared (NIR) and radio observations of the afterglow of GRB 011121. 
We undertake afterglow modeling of this important event and to our delight have found a 
good case for a wind-fed circumburst medium. Thus, the totality of the data — the HST 
optical lightcurves and multi-wavelength (radio, NIR, and optical) data — now support a 
massive star origin for this GRB. 

2. Observations 

Gamma- Rays: GRB 011121 was observed by numerous spacecraft in the Interplanetary 
Network: Ulysses, BeppoSAX (GRBM), HETE-2 (FREGATE), Mars Odyssey (HEND) and 
Konus-Wind. The Tgo duration, as determined from the Ulysses data, was 28 s, placing this 
event in the class of "long bursts" (Figure 1). The peak flux in the 25-100 keV range, over 
0.25 s, was 2.4 x 10~^ ergcm~^s~^, and the fiuence was 2.4 x 10^^ ergcm^. 

Near-Infrared: We observed the afterglow in the near-infrared with the newly-commissioned 
IRIS2 on the Anglo- Austrahan Telescope (AAT), WFIRC on the du Pont 2.5-m telescope 
and the IRCam on the Walter Baade 6.5-m telescope in J and Kg filters. The images were 
dark-subtracted, flat-fielded, sky-subtracted and combined using DIMSUM (Eisenhardt et 
al. 1999) in IRAF. PSF-fitting photometry of the afterglow using DAOPHOT (Stetson 1987) 
was performed relative to point sources in the field. Our multiple calibrations are consistent 
with each other and we estimate the systematic error to be less than 0.05 mag (see Table. 1). 

Radio: We initiated observations of GRB 011121 with the Austraha Telescope Compact 
Array (ATCA) starting on 2001 November 22.58 UT (see Table 2). The data were re- 
duced and imaged using the Multichannel Image Reconstruction, Image Analysis and Display 
(MIRIAD) software package. 
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3. Modeling the Afterglow 

3.1. Dust Extinction 

In Figure 2 we display the optical/NIR spectrum of GRB 011121. The apparent cur- 
vature in the spectrum indicates a large magnitude of dust extinction. In view of this, 
estimating the dust extinction accurately is critical not only for the afterglow modeling but 
also as an important input parameter for the supernova modeling of the HST lightcurves 
(Bloom et al. 2002). 

Prom the IR dust maps (Schlegel, Finkbeiner & Davis 1998) we estimate Ay 1.6 mag. 
However, the IR maps have low angular resolution. Indeed, it appears that the line-of-sight 

to the afterglow passes through the edge of a dust cloud ~ 45 arcmin in extent. Fortunately, 
the availability of both the optical and NIR afterglow data allow us to directly estimate the 
extinction along this line of sight directly. 

We make the reasonable assumption that the optical/NIR afterglow follows the standard 
power-law model, cx t^"z/~^, and we apply the parametric extinction curves of Cardelli, 
Clayton & Mathis (1989) and Fitzpatrick & Massa (1988) along with the interpolation 
suggested by Reichart (2001). Thanks to the abundance of our NIR data, which suffers little 
extinction, we can break the degeneracy between j3 and the magnitude of the extinction, Ay. 

In addition to our own measurements we have included those reported in the literature 
(and noted in Figure 2). Since late-time measurements are increasingly dominated by an 
uncertain mix of the afterglow, the host galaxy and the nearby star B (Bloom et al. 2002) 
we restrict the analysis to data obtained over the first two days. 

Our best fit has an unacceptable = 66 for 48 degrees of freedom, but this is mainly 
dominated by outliers, particularly in the data from the AAT where the seeing blended star 
B with the afterglow in some observations. Inserting an additional 3% error decreases the 
to match the number of degrees of freedom. The additional error term, while ad hoc, is 
reasonable given the variety of telescopes and reduction techniques in our data set. 

Our measured extinction is Ay = 1.16 ± 0.25 mag, distinctly lower than that deduced 
from the dust maps. The type of extinction curve (e.g.. Milky Way, LMC, SMC etc.) is 
unconstrained by these observations. We have not solved for extinction within the host 
galaxy, but the off-center location of the GRB (Bloom et al. 2002) makes it likely that 
the contribution from extinction within the host galaxy is small. Finally, we measure a — 
1.66 ± 0.06 and /3 = 0.76 ± 0.15, without assuming any specific afterglow model. 
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3.2. Afterlow Models 

Armed with a and fi we now consider three afterglow models: (i) isotropic expansion into 
a homogeneous medium (Sari, Piran & Narayan 1998), (ii) isotropic expansion into a wind- 
stratified medium (Chevalier & Li 1999), and (iii) coUimated expansion into a homogeneous 
or wind-stratified medium (Sari, Piran &; Halpern 1999). The models can be distinguished 
by a closure relation, a -|- 6/3 -|- c = 0. These closure relations are due to the dependence of 
both a and /3 on the electron energy distribution index, p, and the values of h and c depend 
on the location of the cooling frequency, i^c, relative to the optical/NIR frequency, vq^ at the 
epoch of the observations. 

As can be seen from Table 3, models with isotropic expansion into an homogeneous 
medium (or, equivalently, a jet which becomes apparent on a timescale longer than the 
epochs of the optical/NIR data used here, tj 2 d) are ruled out by the closure relations at 
more than 2cr significance. Two models produce closure consistent with zero: i^A) A wind 
model with > vq (effective epoch day 1), and p = 2.55 ± 0.08; and (B) A fully developed 
jet at the time of the first optical observation, tj < 0.5 d, with Uc < vq and p = 1.66 ± 0.06. 

The radio measurements, however, do not show any sign of a decay until at least ~ 7 
days after the burst (Figure 3). The rising centimeter-band flux prior to this time indicates 
that the jet break is not at early times, and hence model 5, the jet model, is ruled out. This 
then leaves us with model ^4, the wind model. 

3.3. A Wind Model 

The multi-wavelength data, radio through optical, can only be analyzed by considering 
the evolution of the broad-band synchrotron spectrum. In addition to v^^ p and Ay we need 
to consider the self-absorption frequency, z/q, and the peak frequency, z/^, as well as the 
peak flux, Fi^^rn- These parameters are estimated from the data and can be inverted to yield 
physical quantities, i.e. the energy of the fireball, the density of the ambient medium, and the 
fractions of energy in the electrons, eg, and magnetic field, eb- An example of this approach 
can be found in Berger et al. (2001). The density in the wind model is parameterized by A*, 
which is defined through A — M/AnVw — 5 x lO^^A* g/cm (see Chevaher & Li 1999) where 
Vw is the wind speed and M is the mass loss rate. The normalization of ^4* = 1 applies for 
a typical Wolf-Rayet wind speed of 10^ km/sec and M — 10~^ Mq/yt. 

Given the sparse data we prefer to undertake the model fitting in an evolutionary ap- 
proach rather than performing a blind minimization search. For example, we fix the value 
of p and Av to that determined earlier since the radio data has little bearing on these pa- 
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rameters. Next, we know that Uc> uq, but there are no data in the X-ray band to actually 
constrain the value of Vc- We therefore use i^c ~ 10^^ Hz since this is effectively the lowest 
value the cooling frequency can have in this model. We will, at a later point, revisit this 
issue and examine the consequences of increasing Uc- 

The remaining free parameters'*^^, I'm, and F^^^ are relatively easy to constrain for the 
following reasons. The value of Fi^^m determines the overall scaling in both the optical/NIR 
and radio bands, and is therefore constrained by two sets of data. The value of is 
constrained by the turnover in the radio lightcurves (at t ~ 7 days; see Figure 3), as well 
as the flux density of the optical/NIR lightcurves, since for a given value of Fj/,^, the flux 
density in the optical/NIR band is determined by Um- 

Finally, z/^ is constrained by the spectral slope between the two centimeter bands. The 
comparable fiux between 4.8 and 8.7 GHz suggests that < A.SGHz. An independent 
constraint on is also provided by the equation due to Sari & Esin (2001): 

c . o.oe(i . {^r ii^r ii^r (^) " < 

where 77 = min[(z/c/z/m.)~*'^~^'*'^^, 1] is the fraction of the electron energy radiated away. 

We find F^,™ ^ 3 mJy ^ 10^^ Hz, ^ 2 x 10^^ jjz, z/„ ^ 1.4 GHz and r] = 0.2 
provide an adequate description of the afterglow data. From these parameters we obtain 

~ 0.01 and that inverse Compton cooling is marginally important. Higher values of 
are possible in the inverse Compton-dominated regime and if C << 1, implying that is 
well below the centimeter bands. With these observations, we are unable to constrain such 
a model. 



4. Discussion &; Conclusions 

GRB 011121, a relatively nearby burst (z = 0.36), has shot to fame given what appears 
to be firm identification of an underlying supernova component (Bloom et al. 2002). Here 
we presented early time NIR and comprehensive dual-frequency cm- wave observations of the 
afterglow. Thanks to the NIR data, we have been able to accurately measure the considerable 
Galactic extinction towards the burst. Ay — 1.16 ±0.25 mag, significantly smaller than that 
derived from extrapolations of the IR maps (Schlegel, Finkbeiner & Davis 1998). Our value 



^"^Unless otherwise stated, all time-dependent parameters are evaluated at epoch 1 day e.g. = Vmit = 
Id). 
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of Ay is an important physical parameter in the modehng of the underlying SN component 
(Bloom et al. 2002). 

If indeed long duration gamma-ray events such as GRB 011121 are linked to SNe then the 
progenitors of GRBs are massive stars. Such stars possess strong winds and one expects to 
see a signature of the wind- fed circumburst medium (Chevalier & Li 1999). The optical/NIR 
data alone rule out an isotropic explosion in a constant circumburst medium model. The 
radio data firmly rule out a model in which a jet is fully-developed at t < 0.5 d, but allow 
for a wind-fed circumburst medium. We estimate the mass loss rate, M ^ 10~''f "3 MqJI^^ 
where is the wind speed in units of 10^ km s~^. In the collapsar model (MacFadyen, 
Woosley & Heger 2001), one expects the progenitors of GRBs to be massive stars which 
have lost their hydrogen envelopes, i.e. Wolf-Rayet stars. For such stars, ~ 10^ km s~^. 

Interestingly enough, this mass loss rate is similar to that inferred for the progenitor of 
the Type Ic SN 1998bw, 2.5 x lO-^Moyr"! (Li & Chevalier 1999), based on the analysis of 
the radio hght curves (Kulkarni et al. 1998). This unusual SN is thought to be associated 
with GRB 980425 based on spatial and temporal coincidence (Galama et al. 1998), as well 
as its relativistic outflows (Kulkarni et al. 1998). However, this GRB, if associated with 
SN 1998bw (as we believe), releases at least three orders of magnitude less energy in gamma- 
rays compared to cosmological bursts (Galama et al. 1998) such as GRB 011121. So the 
relation of GRB 980425 to cosmologically located GRBs is unclear. Nonetheless, we make 
the following curious observation: the 7-ray profile (Figure 1) is of similar duration and 
smoothness (with a few spikes superposed) as that of GRB 980425. 

The current data clearly rule out a jet break on the timescale of the optical data, tj ^ 2 
d, and the radio data require tj 7 d. In the formulation of Frail et al. (2001) the opening 
angle of the jet must be wider than 6j > 10 degrees and hence the true energy release is 
larger than 5 x 10^° erg. This lower limit is consistent with the the clustering of energies 
around 5 x 10^'^ erg found by Frail et al. (2001). 

Further improvements to the modeling is possible by including the BeppoSAX measure- 
ment of the X-ray afterglow (Piro et al. 2001); the X-ray fiux will pin down quite well. 
We also note that the radio fluxes given in Table 2 suffer from strong variability (due to 
interstellar scintillation). Here we have used the mean fluxes, and in a later paper we intend 
to report detailed analysis of the scintillation and include the variability as a part of our 
afterglow modeling, in particular as a way to constrain the size of the afterglow region (c.f. 
Frail et al. 1997). 

Thus, at least for one long duration burst the SN-GRB connection and a massive pro- 
genitor origin appears to to have been established. However, the true story may be more 
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complex. The absence of SN components in other GRBs can be explained by appeahng 
to the well known wide diversity in luminosity of Type Ib/c SNe. However, some of the 
intensively observed afterglows are best modeled by expansion into a homogeneous medium. 
There could well be two different classes of progenitors within the class of long-duration 
GRBs (Chevaher & Li 2000). 
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Table 1. NIR observations of the afterglow of GRB 011121. 
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K 


16, 


,537 




0, 


.041 


AAT 


Nov 


22 


.5126 


K 


16, 


,495 


± 


0, 


.035 


AAT 


Nov 


22, 


.6066 


K 


16, 


,605 


± 


0, 


.058 


AAT 


Nov 


22 


.6169 


K 


16, 


,788 


± 


0, 


.042 


AAT 
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Table 1 — Continued 



Date (2001 UT) 


Filter 


Magnitude 


Telescope 


Nov 22.6397 


K 


16.782 ± 0.038 


AAT 


Nov 22.6506 


K 


16.862 ± 0.036 


AAT 


Nov 22.6612 


K 


17.019 ± 0.052 


AAT 


Nov 22.6716 


K 


16.852 ± 0.039 


AAT 


Nov 22.6822 


K 


17.035 ± 0.083 


AAT 


Nov 22.7272 


K 


17.005 ± 0.051 


AAT 


Nov 22.7384 


K 


17.087 ± 0.079 


AAT 


Nov 23.3336 


K 


17.924 ± 0.051 


dP 


Nov 28.7092 


K 


19.346 ± 0.234 


AAT 



Note. — (a) Observations at the du Pont (dP) 2.5-m 
were made by K. Koviak; observations at the AAT were 
made by S.D. Ryder (Nov 22) and K. Gunn (Nov 28); 
observations at the Baade telescope were made by M. 
PhiUips. (b) The following reference stars were used. 
For Kg observations on the AAT we observed UKIRT 
Faint Standards FS 7, 11 and 13 (Hawardcn et al. 2001) 
on 2001 Nov. 28. SJ9113 (Persson, S. E. and Murphy, 
D. C. and Krzeminski, W. and Roth, M. and Rieke, 
M. J. 1998) was observed at the du Pont telescope on 
2001 Nov. 23. We assumed an atmospheric extinction 
coefficient of 0.09 mag/airmass in K for the du Pont 
observations, and that the colour terms were negligible. 
We used the reference stars calibrated by Phillips et al. 
(2001) to cahbrate our J-band observations, (c) The 
AAT measurement of Nov. 28th is contaminated both 
by the host and the nearby star. 
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Table 2. Radio Observations of GRB 011121 made with the Austraha Telescope Compact 

Array. 



Epoch 


1^0 


S±(7 




(UT) 1 


(GHz) 


(A*Jy) 


2001 


Nov 25.20 


4.80 


240 ± 70 


2001 


Nov 28.64 


4.80 


510 ± 38 


2001 


Dec 6.80 


4.80 


350 ± 42 


2001 


Dec 15.80 


4.80 


250 ± 34 


2001 


Dec 22.90 


4.80 


-99 ± 49 


2002 


Jan 23.85 


4.80 


320 ± 38 


2001 


Nov 22.83 


8.70 


210 ± 40 


2001 


Nov 25.08 


8.70 


450 ± 130 


2001 


Nov 28.64 


8.70 


610 ± 39 


2001 


Dec 6.80 


8.70 


220 ± 58 


2001 


Dec 15.80 


8.70 


274 ± 37 


2001 


Dec 22.90 


8.70 


237 ± 46 


2002 


Jan 23.85 


8.70 


-99 ± 47 



Note. — The columns are (left 
to right), UT date of the start 
of each observation, center fre- 
quency, and peak flux density at 
the best fit position of the radio 
transient, with the error given 
as the root mean square noise 
on the image. All observations 
were obtained using the contin- 
uum mode and a 128 MHz band- 
width. Flux calibration was per- 
formed using PKS B1934-638, 
while the phase was monitored 
using PKS B1057-797. 
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Table 3. Afterglow Model Testing 



Model 




(6,c) 


Closure 


P 


ISM 


B 


(-3/2,0) 


1.04 ±0.47 


3.21 ±0.08 


ISM 


R 


(-3/2,1/2) 


2.04 ±0.47 


2.88 ±0.08 


Wind 


B 


(-3/2,-1/2) 


0.04 ±0.47 


2.55 ±0.08 


Wind 


R 


(-3/2, 1/2) 


2.04 ± 0.47 


2.88 ± 0.08 


Jet 


B 


(-2,-1) 


-0.86 ±0.31 


1.66 ±0.06 


Jet 


R 


(-2,0) 


0.14 ±0.31 


1.66 ±0.06 



Note. — Calculation of the closure relations a ± 6/3 ± c 
for a variety of afterglow models. A successful model will 
have a value of zero for the closure relation. The ISM 
and Wind models are for isotropic expansion in an ho- 
mogeneous and wind-stratified medium respectively. The 
Jet model is for coUimated expansion, with the jet break 
time before the first observations were made. The rela- 
tions are dependent on the assumed location of the cooling 
frequency, Uc relative to the optical and NIR bands, vq'- 
the case Vc > vo is denoted by "B" (for blueward) and 
Vc < vq by "R" (for redward). p is the electron energy 
power law index. 
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Fig. 1. — Time history of GRB 011121 in the 25-150 keV energy range, as observed by 
Ulysses. The dashed hne gives the background rate. Zero on the time axis corresponds to 
an Earth-crossing time of 67630.899 s. 
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Fig. 2. — The optical and NIR spectral flux distribution of the afterglow of GRB 011121 
at 0.5 d, based on measurements presented here (circles) and from the literature (triangles; 
Olsen et al. 2001; Stanek & Wyrzykowski 2001). Measurements taken within 0.5 ± 0.1 d 
have been transformed to 0.5 d using the best fit model. The sohd lines indicates our best 
fit to the data, using a power-law model plus foreground extinction. The dashed line is the 
intrinsic spectrum of the afterglow. 
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Fig. 3. — The radio light curve of the afterglow of GRB 011121. The solid line is our wind 
model, the thin line is a representative jet model, which is clearly excluded by the data. The 
radio data exhibit strong modulation due to interstellar scintillation and, as a result, deviate 
from our model by more than la. This will be addressed in a future paper. 
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